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ABSTRACT 

We study transverse fast magnetohydrodynamic waves in a system of two 
coronal loops modeled as smoothed, dense plasma cylinders in a uniform magnetic 
field. The collective oscillatory properties of the system due to the interaction 
between the individual loops are investigated from two points of view. Firstly, 
the frequency and spatial structure of the normal modes are studied. The system 
supports four trapped normal modes in which the loops move rigidly in the trans- 
verse direction. The direction of the motions is either parallel or perpendicular 
to the plane containing the axes of the loops. Two of these modes correspond to 
oscillations of the loops in phase, while in the other two they move in antiphase. 
Thus, these solutions are the generalization of the kink mode of a single cylinder 
to the double cylinder case. Secondly, we analyze the time-dependent problem 
of the excitation of the pair of tubes. We find that depending on the shape 
and location of the initial disturbance, different normal modes can be excited. 
The frequencies of normal modes are accurately recovered from the numerical 
simulations. In some cases, because of the simultaneous excitation of several 
eigenmodes, the system shows beating. 

Subject headings: Sun: corona-magnetohydrodynamics (MHD)-waves 



1. Introduction 

Transverse coronal loop oscillations have been studied in recent years after being ob- 
served for the fi rst time by the Transition Reg i on and Coronal Exp l orer (TRACE) in 199 8 



see for example lAschwanden et al.lll999l . |2002| ; ISchrijver et al.ll2002l ; IVerwichte et al.ll2004j ). 



These oscillations were initiated shortly after a solar flare that disturbed the loops. Much 
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before the TRACE observations , the t heory of loop oscillations was developed (jSpruitlll981 
Edwin and Roberts! Il983l ; ICallyi Il986l ) and the different kinds of oscillations were studied. 
The observed transverse motions have been interpreted in terms of the excitation of the fast 
magnetohydrodynamic (MHD) fundamental kink mode. 

Most analytical studies about transverse loop oscillations have only considered the prop- 
erties of individual loops, but in many cas es loops belong to com plex active regions where 
they are usually not isolated. For example, IVerwichte et al.l (12004 ) reported complex trans- 
verse motions of loop s in a post-flare arcade. In particular, loops D and E (see Fig. 1 of 
Verwichte et al.ll2004j ) show bouncing displacements with oscillations in phase and antiphase 
that repeat in time. The same behavior of the movements in a loop bundle c an be observed 



in th e event of March 23, 2000 of the compact flare recorded by TRACE (see lSchrijver et al. 
2002h . Additionally, antiphase oscillations of adj acent loops have also been reported in 



Schrijver and Brown I (120001 ) ; ISchrijver et al.l (120021 ) . These observations suggest that there 
are interactions between neighboring loops and that the dynamics of loop systems is not 
simply the sum of the dynamics of the individual loops. 

On the other hand, it is curren tly debated whethe r active region c oronal loops ar e 
monolithic (lAschwanden et al.l 120051 ) or multi-stranded (iKlimchukl 120061 ; iDeForestl 120071 ) . 
The strands are considered as mini-loops for which the heating and plasma properties are 
approximately uniform in the transverse direction. In the multi-stranded model it is sug- 
gested that loop s are formed by b undles of several ten s or se veral hundreds of physically 
related strands (IKlimchukl 120061 ) . lLopez Fuentes et al.l (120061 ) suggest that these strands 



wrap around each other in complicated ways due to the random motion of the foot points 
in the solar surface. These models explain the constant width and symmetry of the loops as 
observed with current X-ray and EUV telescopes. 

From the observations, it is thus necessary to study not only individual loops but also 
how several loops or strands can oscillate as a whole, since their joint dynamics can be dif- 
ferent from those of a s ingle loop. Little work has been done on composite structures so far. 
Berton and Heyvaertsl (119871 ) studied the magnetoh ydrodynamic normal mo d es of a periodic 



magne tic medium, while other authors, for example iBogdan and Fox! (Il99ll ); iKeppens et al. 



(119941 ). analyzed the scattering and absorptio n of ac o ustic waves by bundles of ma gnetic 
flux tubes with sunspot properties. iMurawskil (119931 ); iMurawski and Roberts! (119941 ) stud- 
ied numerically the propagation of fast waves in two slabs unbounded in the longitudinal 
direction. On the other hand, in iDiaz et al.l (120051 ) the oscillations of the prominence thread 
structure were investigated. These authors found that in a system of equal fibrils the only 
non-leaky mode is the sy mmetric one , with all fibrils oscillating in spatial phase with the 
same frequency. Finally, iLuna et al.l (120061 ) found that in a system of two coronal slabs, 
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the symmetric and antisymmetric modes can be trapped and that an initial disturbance can 
excite these modes, which are readily detectable after a brief transient phase. If the funda- 
mental symmetric mode and the antisymmetric first harmonic are excited at the same time, 
a beating phenomenon takes place. In such a case, the loops interchange energy periodically. 
In any case, all these authors found that a system of several loops behave differently from 
an individual loop. 

Here we consider a more complex system than those studied in previous works. Our 
model consists of two parallel cylinders, without gravity and curvature. This model allows us 
to study the interaction between loops and the collective behavior of the system. We study 
the normal modes and also solve the time-dependent problem of the excitation of transverse 
coronal loop oscillations. We concentrate on a planar pulse excitation and compare the 
results of the simulations with the eigenmodes of the configuration. 

This paper is organized as follows. In £J2] the loop model is presented. In §|3] the normal 
modes are calculated and the frequencies and spatial distribution of the eigenfunctions are 
studied. The time-dependent problem is considered in §U where the velocity and pressure 
field distribution are analyzed for different incidence angles of the initial perturbation. In 
§5] the loop motions are studied and the beating is analyzed. Finally, in §|6] the results are 
summarized and the main conclusions are drawn. 



2. Equilibrium configuration and basic equations 

The simplest way to investigate the int eraction of a set of loops is to consider a pair of 



loops in slab geometry. In lLuna et al.l (120061 ) this model was studied in detail using the ideal 
MHD equations and the zero-/? plasma limit. Here a more realistic model is considered. The 
equilibrium configuration consists of a system of two parallel homogeneous straight cylinders 
of radius a, length L, and separation between centers d (see Fig. [[]). We assume the following 
equilibrium plasma density profile: 



Po(x,y) 



p e , if r± > a and r 2 > a, 
Pi, if T\ < a or r 2 < a, 



where x, y are the Cartesian coordinates and r\ and r 2 , defined as = (x + d/2) 2 + y 2 
and r\ = (x — d/2) 2 + y 2 , are the distances from the point (x,y) to the centers of the left 
and right loops, respectively. In the previous expression p e and p; are the densities in the 
external medium or corona and the loop (p ; > p c ), respectively. Hereafter, we use a density 
contrast p;/p e = 10. 




Fig. 1. — Sketch of the model, that consists of two straight cylinders immersed in the coronal 
medium. Hereafter the total pressure and the velocity fields are plotted in the xy-plane, shown as 
a white slice. 
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The loop centers lie on the x-axis at x = d/2 for the right loop and x = —d/2 for the left 
loop. The configuration is symmetric with respect to the yz-plane and the z-axis is parallel 
to the axes of the cylinders. The tubes and the environment are permeated by a uniform 
magnetic field along the ^-direction (B = B e z ). The Alfven speed, va = B / y/JIp, takes 
the value t>Ai inside the loop and v^e in the surrounding corona (v^i < t>A e )- 

Linear perturbations about this equilibrium for a perfectly conducting fluid in the zero-/? 
limit can be readily described using the ideal MHD equations in Cartesian coordinates. The 
velocity is denoted by v = (v x , v y , 0) and B = (B x , B y , B z ) is the magnetic field perturbation. 
We have assumed a z-dependence of the perturbations of the form e~ lkzZ . In this model we 
consider the photosphere as two infinitely dense planes located at z — ±L/2. The loop feet 
are anchored in these planes and so the fluid velocity is zero at these positions (this is the 
so-called line-tying effect). This condition produces a quantization of the z-component of 
the wave- vector to k z = nir/L. Hereafter we concentrate on the fundamental mode, with 
n = 1. The total pressure perturbation is 

Pt = —B z , (1) 

and coincides with the magnetic pressure perturbation in the zero-/? limit. 



3. Normal modes 



Analytical solutions to the eigenvalue problem of the previous model (assuming a tempo- 
ral dependence of the form e luJt ) are very diffi cult to derive due to the ge ometry of the system. 
The methods used for a single cylinder (see Edwin and Roberts! Il983l ) cannot be applied to 
the st udy of two tubes. One way t o solve the problem is to us e scattering theory, see for ex - 
ample lEdwin and Roberts! (119831). iBoedan and Knolkerl (Il99ll ). iBogdan and Zweibel Jl985h . 
Bogdan and Fox! (Il99ll ) and iKeppens et al.l (119941 ) . Another way is to solve the eigenvalue 
problem given by the ideal MHD equations numeric ally. We hay e used this approach and we 
have done the computations with the PDE2D code (jSewell!l2005l ). We have used bicylindrical 
orthogonal coordinates, defined by the transformation 



x 



d/2 sinhf 
cosh v — cos u 



d/2 sinii 
cosh v — cos u 



(2) 



where < u < 2tt and — oo < v < oo. The loop boundaries are coordinate surfaces 
at v = iarcsinh where the positive and negative signs correspond to the right and 
left tubes, respectively. We impose the restriction that the solutions tend to zero at large 
distances from the cylinders, i.e. we seek trapped mode solutions. 
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We find four collective fundamental trapped modes (see Fig. [2]). There are other har- 
monics but we concentrate on the fundamental kink-like modes because they produce the 
largest transverse displacement of the loops axes. The velocity field is more or less uniform 
in the interior of the loops, and so they move basically as a solid body, while the external 
velocity field has a more complex structure. The four velocity field solutions have a well 
defined symmetry with respect to the y-axis. In Figure [2^, we see that the velocity field 
inside the tubes lies in the ^-direction and is symmetric with respect to the y-axis. We 
call this mode S x , where S refers to the symmetry of the velocity field and the subscript x 
refers to the direction of the velocity inside the tube. The same nomenclature is used for 
the other modes. In Figure [2b the velocity inside the cylinders is mainly in the x-direction 
and antisymmetric with respect to the y-axis, so we call this mode A x . Similarly, in Figure 
[2b the velocity lies in the y-direction and is symmetric with respect to the y-axis, while it 
is antisymmetric in Figure [2H. Hence, we call these modes S y and A y , respectively. The 
pressure field of the A x and S y modes is symmetric with respect to the y-axis, while that of 
the S x and A y modes is antisymmetric. 
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Fig. 2. — Total pressure perturbation (color field) and velocity field (arrows) of the fast 
four collective normal modes (plotted in the xy-plane, see Fig. [1]). The modes are labeled 
as (a) S x , with the loops moving in phase in the x-direction; (b) A x , the tubes move in the 
x-direction but in antiphase; (c) S y , the tubes move in the ^-direction in phase; and finally 
(d) A y , the loops move in antiphase in the ^/-direction. Here, the loop radii are a = 0.1L 
and the distance between centers is d — 6a. 
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The frequencies of oscillation of these four modes as a function of the loop separation, 
d, are displayed in Figure [3) For large separations between the tubes, the modes tend to the 
kink mode of an individual loop (see dotted line). On the other hand, for smaller separations, 
they split in four branches a ssociated to the fo ur o scillatory modes d escribed before. The 



splitting effect was noticed in iDiaz et al.l (120051 ) and iLuna et al.l (120061 ) in a configuration of 



several slabs. The frequency difference between the modes increases when the interaction 
between the loops becomes stronger, i.e. when the distance between them is small. When 
the loops are very close (d ~ 2a), the frequencies of the S x and A y modes tend to the value 
L) = 3.33/rAi, which is similar to the internal cut-off frequency, uj c \ = k z VAi = 3.14/rAi (the 
difference is only around 6%). Here tm is the Alfven transit time, defined as ta; = L/vj^. 
On the other hand, in this limit, the S y and A x frequencies are quite large in comparison to 
the kink mode frequency. 

It is interesting to note that when both tubes move symmetrically in the ^-direction, 
i.e. in the S x mode, the fluid between follows the loops motion (see Fig. Wfr)- On the other 
hand, when the loops oscillate antisymmetrically, i.e. in the A x mode, the intermediate fluid 
is compressed and rarefied (see Fig. [2b), producing a more forced motion than that of the 
symmetric mode. This is the reason for the S x (A x ) mode having a smaller (larger) frequency 
than that of the individual loop. For the modes polarized in the ^/-direction the behavior 
is somehow similar, although in this case the antisymmetric mode (see Fig. [2H) has a lower 
frequency than the symmetric mode (see Fig. [2fc). When one of the loops moves upwards 
the surrounding fluid near the other loop moves downwards. This helps to push the other 
loop in this direction and produces the antisymmetric motion. The situation is different for 
the S y mode, for which the direction of motion of the surrounding fluid is opposite to that 
of the other tube. This explains why the frequency of the A y solution is smaller than that 
of the S y mode. 
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Fig. 3. — Frequency, u>, as a function of the separation between cylinders, d, for a density 
enhancement pi/p — 10 and loop radius a = 0.1L. The lower horizontal thick dotted line 
is the internal cut-off frequency, uj c \ = 3.14/tai- The thin dotted line is the kink frequency 
(u = 4.12/rAi) of an individual loop. The calculated frequencies from the time-dependent 
results in §4.11 are also plotted as triangles (S x ), squares (A y ), asterisks (S y ), and diamonds 
(A x ). 
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4. Time-dependent analysis: numerical simulations 

The initial perturbation that we have used when solving numerically the ideal MHD 
equations is a planar pulse in the velocity field of the form 

y = Ke^ K < T -^^\ (3) 

i.e. a Gaussian profile (of width A = 0.2L centered at r = (d/2,0, 0)) and direction of 
propagation along K = — (cos a, sin a, 0), a being the angle between the wavevector and the 
x-axis. Here K also defines the initial polarization of v, which is perpendicular to the planar 
pulse. The initial value of the magnetic field perturbation is zero, and thus the same applies 
to the total pressure perturbation. In the simulations a spatial domain of size 30a x 30a is 
used and the boundaries are far from the loops. These boundaries are open, which ensures 
that the numerical reflections are negligible. 

In Figures HI and [6] three examples of the time evolution are shown for a = 90°, 0° 
and 45°, respectively, and for a fixed distance between loops d = 6a, identical to the one 
used in Figure [2] (see the time evolution in Movie 1, Movie 2, and Movie 3). These three 
cases illustrate the time evolution of the system after a perturbation, which consists of two 
regimes: the transient and the stationary phases. The stationary phase is characterized by 
oscillations in one or several fundamental trapped normal modes (see $3]). On the other 
hand, in the transient phase there are leaky modes and internal reflections and refractions. 
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Fig. 4.— Time-evolution of the velocity field (arrows) and total pressure field (colored 
contours), for a separation between loops d = 6a and an initial pulse with an angle a = 90°. 
The two circles mark the positions of the loops at t = 0. The panels show different evolution 
times. In (a) the initial condition over the velocity field is represented. In (b) the velocity 
and pressure field shortly after the initial disturbance, that is, during the transient phase, 
are shown. Both tubes are excited at the same time. In panel (b) the tubes are in the 
transient phase. In panels (c) and (d) the system oscillates in the stationary phase with the 
S y normal mode. This time evolution is also available as an mpeg animation in Movie 1. 
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In Figure H] (see Movie 1) the time evolution for the a = 90° initial disturbance is shown, 
for which, the pulse front lies along the x-axis and excites the v y component. The loops are 
perturbed at the same time (as can be appreciated in Fig. Hb) and as a consequence they 
oscillate symmetrically. In Figure Hb the system is in the transient phase, characterized 
by internal reflections related with the emission of leaky modes. The external medium has 
not relaxed yet. Finally, the system reaches the stationary phase (see Figs. [2b and HH) and 
oscillates with the S y trapped mode (compare the velocity field and the pressure distribution 
with Fig. [2b). 



13 




Fig. 5. — Same as Figure H] for an initial pulse with an angle a = 0°. Here the stationary 
phase is governed by a superposition of the S x and A x normal modes. The whole time 
evolution is presented in Movie 2. 
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In Figure [5] (and Movie 2), the time evolution for the a = 0° initial disturbance is shown. 
Now the pulse is centered on the right loop (see Fig. [5^) and excites the v x component. In 
Figure Eb, the pulse reaches the left tube and passes through it, the system still being 
in the transient phase. On the other hand, in Figures \5jp and EH the system oscillates in 
the stationary phase. It is interesting to note that this particular initial disturbance does 
not excite the left loop; neither at t = nor during the transient phase. Nevertheless, 
the oscillatory amplitude in the left loop grows with time in the stationary phase, while 
the amplitude in the right loop decreases in the time interval shown in Figures \5jp and EH 
(see also Movie 2). Then, it is clear that the left tube acquires its movement through the 
interaction with the right loop, i.e. by a transfer of energy from the right loop to the left 
loop. This process is reversed and repeated periodically: once the left loop has gained most 
of the energy retained by the loops system, so that the right loop is almost at rest, the 
left tube starts giving away its energy to the right cylinder, and so on. This is simply a 
beating phenomenon, that can be explained in terms of the normal modes excited in this 
numerical simulation. In fact, the initial disturbance excites the S x and A x modes with the 
same amplitude and for this reason the excitation is initially maximum on the right tube 
and zero on the left tube. A more detailed discussion about this issue is given in §5J 
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Fig. 6. — Same as Figure H] for an initial pulse with an angle a = 45°. Here the stationary 
phase is governed by a superposition of the S x , A x , S y , and A y normal modes. The whole 
time evolution is presented in Movie 3. Movie 4 contains the time evolution for much larger 
times. 
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Finally, we discuss the results for an excitation with a = 45°. This simulation is the most 
complex and general of all (see Movie 3). As we can see in Figure 16^, now both components 
of the velocity are excited. In Figure Eb the initial pulse reaches the left tube and passes 
through it, but only leaky modes are excited. In Figures and EH the system oscillates in 
the stationary phase, which is a combination of the four modes S x , A x , S y and A y . As in the 
previous case, there is beating but now it is present in both the x- and y- velocity components. 
Like for the previous simulation, the left loop is almost still until the stationary phase (see 
also dotted curves in Figs. [7^ and [7H) despite that in this simulation the pulse directly hits 
the left loop without the obstacle of the right loop. In §5] details about the behavior of the 
system are given. 

Once we know the general features of the excitation of the two cylinders we can perform 
a parametric study of the effect of the distance between the loops and also the angle of 
excitation on the loops motion. 

4.1. Effect of the distance between loops 

We generate an initial disturbance with an angle of 45° for different distances d and 
measure the velocity in the loops as a function of time. From this information we can 
extract the frequencies of oscillation. As we have seen, since the velocity field inside the 
loops is more or less uniform (see Fig. [2]), it is enough to measure the velocity at the center 
of the loops to describe their global motion. The reason for choosing the initial disturbance 
with a = 45° is that it excites the four normal modes, so that with a single simulation we 
can measure their frequencies. 
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Fig. 7. — (a) x-component and (d) y-component of the velocity at the center of the right 
(solid line) and left (dotted line) loops for the numerical simulation of Figure [6] (i.e. with 
an initial incidence angle a = 45°). With the method explained in §4.11 the normal mode 
velocities are extracted and plotted in (b) for the S x (solid line) and the A x (three-dot-dashed 
line) modes and in (e) for the S y (dashed line) and A y (dot-dashed line) modes. The corre- 
sponding power spectra are plotted with the same line styles in (c) and (f). Power maxima 
allow us to determine the frequency of the normal modes from the numerical simulation. 
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Fig. 8. — Same as Figure [7] for an initial incidence angle a = 70°. 
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In Figures [7^ and [7U the x- and ^-components of the velocity at the center of each loop 
are plotted. In these figures we see that, after a very brief transient characterized by short- 
period oscillations, the system oscillates with the sum of normal modes. The frequencies 
of the modes are quite similar, and it is difficult to resolve them. Although the frequencies 
of these modes are present in the time- dependent signal, this information cannot be easily 
extracted from the data because in these simulations the maximum evolution time (which 
is determined by the numerical damping) is T = 6t^. With this maximum time we have a 
frequency resolution 2/T ~ 0.3/r Ai , but, as evidenced by Figure El the difference in frequency 
between the eigenmodes is typically less than 0.1/tai so we have not enough frequency 
resolution. For this reason we extract the frequencies with another method considering 
that the velocity field is the addition of normal modes with symmetric and antisymmetric 
spatial functions with respect to the y-axis. We measure the velocity in the loop centers 
(x = —d/2,y = 0) and (x = d/2,y = 0), i.e. two symmetric points with respect to x — 0. 
Then, the sum of both measured velocities in these points is twice the symmetric velocity. 
Dividing this velocity by two we obtain the v x of the S x mode and the v y of the S y mode 
in these points. On the other hand, the subtraction of the measured velocities is twice the 
antisymmetric velocity. Similarly, dividing this velocity by two we obtain the v x of the A x 
mode and the v y of the A y mode. The obtained mode velocities are plotted in Figures [7b 
and [Tfe. Next, we compute a periodogram of these signals (plotted in Figs. and[7J), from 
which the frequencies of the collective modes are determined. The periodogram is preferred 
over the FFT as it allows to precisely identify these frequencies. 

The above procedure has been applied to numerical simulations for different separations 
between loops and the frequencies of the four fundamental eigenmodes have been obtained. 
The calculated frequencies have been superimposed to the normal mode values in Figure [3] 
using symbols. A good agreement between the normal mode calculations and the time- 
dependent results can be appreciated. 

4.2. Effect of the incidence angle 

We next study the evolution of the system for different incidence angles, a, of the planar 
pulse and a fixed distance between loops (d = 6a). Some examples of the time evolution 
have already been discussed and shown in Figures HI El and El The mode excitation depends 
on the width, A, the incidence angle, a, and the position, r , of the initial disturbance, 
but here we only consider the dependence on the incidence angle. The angles considered in 
our simulations vary from a = 0° to 90° with steps of 5°. Using the method of §4.11 it is 
also possible to extract the amplitude of each normal mode, given by the amplitude of the 
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sinusoidal oscillations in the stationary phase. Two examples of the extraction method are 
plotted in Figure [7] for a = 45° and Figure [8] for a = 70°. 




Fig. 9. — Amplitudes of the four normal modes as a function of the incidence angle a. The 
separation between loops is d — 6a. 
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In Figure [9] the amplitude of the four collective modes is plotted as a function of the 
incidence angle. The behavior of the amplitude can be divided in two parts, namely for 
0° < a < 50° and for 50° < a < 90°. In the first interval the amplitudes of the S x and A x 
modes are more or less equal (see Figs. [7b and [Tb as an example) and can be approximated 
by 0.3 cos a. The same occurs for the amplitudes of the S y and A y modes, which vary roughly 
as 0.29 sin a. In the second interval these amplitudes can be quite different (see Figs. [Hb 
and [8k as an example) and the S x , A x , and A y amplitudes go to zero at a = 90°. On 
the other hand, the S y amplitude increases and reaches its maximum value at a = 90°. 
Furthermore, for a = 0° the amplitudes of the S x and A x modes have a maximum around 
0.3 while the amplitudes of S y and A y modes are zero. This is because for a = 0° the initial 
disturbance drives the x-component of the velocity and so only the S x and A x modes are 
excited. Similarly, for the perturbation with a = 90° only the S y and A y modes can be 
excited, although the shape of our initial perturbation prevents the A y mode from being 
driven and so the S y mode reaches the largest amplitude of all modes. On the other hand, 
the excitation of the antisymmetric modes requires the initial disturbance to hit the right 
and left loops at different times. For this reason, the amplitudes of these modes decrease 
with a. In fact, when a = 90° this time difference is zero since both loops are excited at the 
same time and the amplitude of the A x and A y vanishes. Finally, it is interesting to note 
that for a = 45° the four modes are excited with almost the same amplitude. 



5. Study of the loops motions: beating 

As we have shown in the previous section, loop motions can be very complex. This 
is even more clear in Movie 4, in which the time-evolution for a simulation with identical 
parameters to those used in Figure El but for much larger times is represented. In §4] we 
mentioned that the initial disturbance excites the right loop but does not perturb the left 
loop. After a short time the left tube starts to oscillate due to the interaction with the 
right one. At this stage, the right loop oscillates with the velocity polarization of the initial 
pulse, whereas the left tube oscillates in a direction perpendicular to that of the initial 
disturbance. The reason for the complexity of the loops motions is that their oscillations 
are a superposition of four normal modes with different velocity polarizations, parities, and 
frequencies. 
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Fig. 10. — Temporal variation of the velocity components (a) v x and (b) v y at the center 
of the right loop (solid line) and left loop (dashed line). These results correspond to the 
simulation shown in Figure [6] and illustrate the beating of the pair of loops. Damping caused 
by numerical dissipation causes a slight decrease of the amplitude during the numerical 
simulation. The time-evolution is also available as an mpeg animation in Movie 4. 



We next analyze this case in detail. The x- and ^/-components of the velocity at the 
center of the loops are represented in Figures [TUh and [TUb . respectively. There is a clear 
beating, characterized by the periodic interchange of the x- and y- components of the velocity 
between the loops. The two velocity components are modulated in such a way that v x reaches 
its maximum value in the left tube and becomes zero in the right tube at the same time 
(around t ~ 40rAi). This process is reversed at t ~ SOtaz and repeats periodically. 

The loops motions can be studied theoretically. In the stationary phase, during which 
the system oscillates in the normal modes S x , A x , S y , and A y , the velocity field components 
are 

v x (x,y,t) = C s x (x,y)cos(u s x t + ( p s x ) + C£(x,y) cos (u£t + <f>£) , (4) 
v y (x, y, t) = C y s (x, y) cos (w*t + <p s y ) + C^{x, y) cos (uft + <ft) . (5) 

The S and A superscripts refer to the symmetric and antisymmetric modes, respectively. The 
functions Cj, C£, Cy, and Cy represent the spatial distribution of the four normal modes 
(see Fig. |2J) and their amplitude accounts for the energy deposited by the initial disturbance 
in each of them. The normal mode frequencies are represented by their frequencies, u, while 
(J) are their initial phases. 

Let us turn our attention to the results in Figure [7J In the loops centers the symmetric 
and antisymmetric modes have a very similar amplitude (see also Fig. [10] for a = 45°), 
which means that Cf (d/2, 0) = C x (d/2,0). Then, taking into account the parity of Cf and 
C x about x = 0, we have C|(— d/2, 0) = — C£(— d/2,0). Inserting these expressions into 
equations PJ and ([SD evaluated at the loop centers we obtain 

v„ ght («) = (cUos(^±^^ t) cos(^^^^ 

v*(t) = - (C. Bin(^^ t) Bi„(^i t), C v t) AA<±^ ()) .(7) 

where v right and vi eft are the velocity of the right and left loop, respectively. We have defined 
C x = 2Cf (d/2, 0) and C y = 2Cj(rf/2, 0) and have assumed = 4> x = <fi y = <fi y = because 
the initial disturbance is over the right loop. The beating curves shown in Figure [TD] are 
accurately described by these equations. 

These formulae contain products of two harmonic functions. Then, the temporal evo- 
lution during the stationary phase is governed by four periods: the two oscillatory periods, 
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giving the mean periods of the time signal; and two beating periods, 
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giving the periods of the envelop of the time signal. To apply these expressions to the 
numerical simulation of Figure [7] we insert the values of cuf , Uy, and ujy for d = 6a 
into equations (l8l)-(fTTT). Then we obtain T x = 1.52rAi, T y = 1.52rAi, Tta = 159.96tai, and 
T hy = 479.88rAi- The two oscillating periods are equal because the frequency distribution 
is approximately symmetric around the central value (the kink frequency of an individual 
loop) for sufficiently large distances (see Fig. The two beating periods derived from the 
numerical simulations match very well these values because Figure fTUl gives T^ x /4 ~ 40rAi 
and T h y/4 ~ 120r Ai . 

The Ti / 2 phase difference between v right and v left (see Figs. [7^ and [7H) is due to the 
fact that our system of two loops basically behaves as a pair of driven-forced oscillators. 
Considering v x , the left loop has initially a 7r/2 delay with respect to the right loop because 
it behaves as a driven oscillator and the left one like a forced oscillator. After half beating 
period, T^ x /2, the roles are exchanged and left loop becomes the driver and right one the 
forced oscillator. The ^-co mponents of y right a nd v left exhibit the same behavior (see Fig. (TU). 



This was already shown by iLuna et al.l (120061 ) in the case of two slabs. 



As we have seen, the polarization of the oscillations changes with time (see Movie 4 for 
an example). In the beating range, we can see this from the equations by calculating the 
scalar product of the velocity at the loop centers, 

v r ight • vi eft = -C 2 X sin (2(u£ ~ sin (2(u£ + 

-C 2 y sin (2^ - u s y )?) sin (2^ + u s y )t) . (12) 

This product gives the relative polarization of the loop oscillations and we see that is zero 
at t — and approximately zero for sufficiently small times. Thus, the left loop does not 
oscillate initially and it starts to oscillate perpendicularly to the right loop during the first 
oscillations. This feature is shown in Figure [6] and Movie 3 and Movie 4. 

Similar beating features are recovered for incidence angles of the initial disturbance in 
the range 0° < a < 50° (what we call the beating range). The cause of this behavior is 
explained by Figure [9j for these values of a a similar amount of energy is deposited in the 
S x and A x modes, so the beating of the v x component is possible. Obviously, an analogous 
argument applies to v y . This is not the case for 50° < a < 90° for which the symmetric 
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and antisymmetric modes receive different amounts of energy from the initial excitation and 
then their relative amplitude is different (see Fig. [8] for an example). Simulations for angles 
a > 50° do not clearly exhibit beating and the trajectories of the loops are much more 
complex than those in the beating range. 



6. Discussion and conclusions 

In this work, we have investigated the transverse oscillations of a system of two coronal 
loops. We have considered the low-/?, ideal MHD equations and have studied both the normal 
modes of this configuration and the time-dependent problem. The results of this work can 
be summarized as follows: 

• The system has four fundamental normal modes, somehow similar to the kink mode of 
a single cylinder. These modes are collective, i.e. the system oscillates with a unique 
frequency, different for each mode. When arranged in increasing frequency the modes, 
are S x , A y , S y , and A x , where S(A) stand for symmetric (antisymmetric) velocity 
oscillations with respect to the plane in the middle of the two loops and x (y) stands 
for the polarization of motions. These modes produce transverse motions of the tubes, 
so they are kink- like modes. 

• We have studied the eigenfrequencies as a function of the separation of loops. For 
large distances between cylinders, they behave as a two independent loops, i.e. the 
frequency tends to the individual kink mode frequency. When the distance decreases 
the frequency splits in four branches, two of which correspond to the S x and the A y 
modes and are below the frequency of the individual tube, and the other two are 
related to the S y and A x modes and lie above the kink frequency of a single tube. 
Roughly speaking, there is a certain parallelism between our system of two loops and 
a mechanical system of two coupled oscillators with / degrees of freedom, which has 
/ x n collective normal modes. This parallelism is possible because a slab or a cylinder 
oscillating with the kink mode moves more or less like a solid body. The number of 
translational degrees of freedom one for an individual slab (/ = 1) and two for an 
individual loop. Then, the parallelism with the mechan ical system tells us that in a 



two slab system there are two collective normal modes (ILuna et al.l 120061 ). while in a 
two cylinder system there are four. 

For small distances between the loops, the frequency of the S x and A y modes is quite 
similar and tends to the intern al cut-off frequen cy. This is different to the behavior in a 



configuration of two slabs (see ILuna et al.ll2006l ) where, for small distances between the 
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slabs, the system behaves as an individual loop of double width. On the other hand, for 
the two cylinders the frequency is much lower than that of a loop with double radius. 

• We have also studied the temporal evolution of the system after an initial planar pulse. 
We have shown that, depending on the incidence angle, the system oscillates with a 
combination of several normal modes. The frequencies of oscillation calculated from 
the numerical simulations agree very well with the normal mode eigenfrequencies. 

• In the beating range (0° < a < 50°), the system beats in the x- and ^/-components 
of the velocity and the left and right loops are 7r/2 out of phase for each velocity 
component. They behave as a pair of driven-forced oscillators, with one loop giving 
energy to the other and forcing its transverse oscillations. The role of the two loops 
is interchanged every half beating period. On the other hand, for perturbations with 
a > 50° the loops motions are much more complex than those in the beating range. 
The phase lag cannot be clearly appreciated and it strongly depends on the incidence 
angle of the initial pulse. 

From this work, we conclude that a loop system clearly shows a collective behavior, its 
fundamental normal modes being quite different from the kink mode of a single loop. These 
collective normal modes are not a combi nation of individua l loop modes . This suggests 



that t h e observed oscil l ations reported in lAschwanden et al.l (119991 . |2002j) ; ISchrijver et al. 



(120021 ); IVerwichte et al.l (12004 ) are in fact caused by one or a superposition of s ome c ollective 
modes. Moreover, the antiphase movements reported by iNakariakov et al.l (119991 ) can be 
easily explained using our model. The same applies to the bounce movement of loops D 
and E studied in IVerwichte et al.l (120041 ) . These motions can be interpreted by assuming 
that there is beating between the loops produced by the simultaneous excitation of the 
fundamental S x and A x modes. 

It should be noted that the observations indicate a very rapid damping of transverse 
oscillations, such that in a few periods the amplitude of oscillation of the loops is almost 
zero. This fast attenuation may hide the beating produced by the simultaneous excitation 
of several normal modes of the system. However, in some situations, for example, for small 
loop separations and high density contrast loops, the beating periods decrease. Then, under 
such conditions the beating could be detectable in the observation interval. In any case, the 
beating is just one particular collective behavior, and there is always interaction between the 
individual loops in short time scales (typically of the order of 2d/vxe)- The consequences of 
this interaction are the collective normal modes of the system. The presence of the normal 
modes could be also clear from a frequency analysis. Unfortunately, due to the temporal 
resolution, these observations do not allow us to perform such analysis, but the frequency 
extraction method derived in §4.11 is suitable to be applied to the observations. 
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Finally, in order to have more realistic models additional effects need to be included. 
In this work, we have studied two loops with exactly the same density and radii, so the 
next step is to analyze the behavior of a system of n loops with different properties. This 
study could also be extended to understand the possible effect of internal structure (multi- 
stranded models and small filling factors) on the oscillating loops by considering a set of very 
thin tubes with different densities and radii. We expect that the dynamical behavior and 
frequencies of multi-str anded loops di f fer fro m those of the monolithic models. Preliminary 



work has been done by lArregui et al.l (120071 ) who have studied the effects on the dynamics 
of the possibly unresolved internal structure of a coronal loop composed of two very close, 
parallel, identical coronal slabs in Cartesian geometry. 
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